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ABSTRACT 



"^" Context. Stars whose mass is an order of magnitude greater than the Sun play a prominent role in the evolution of galaxies, exploding 

^-H as supernovae, triggering bursts of star formation and spreading heavy elements about their host galaxies. A fundamental aspect of 

star formation is the creation of an outflow. The fast outflow emerging from a region associated with massive star formation in the 
r— I Orion Molecular Cloud 1 (OMC-1), located behind the Orion Nebula, appears to have been set in motion by an explosive event. 

-^F Aims. Here we study the structure and dynamics of outflows in OMC-1. We combine radial velocity and proper motion data for 

Onear-IR emission of molecular hydrogen to obtain the first 3-dimensional (3D) structure of the OMC-1 outflow. Our work illustrates 
^ a new diagnostic tool for studies of star formation that will be exploited in the near future with the advent of high spatial resolution 

r^ spectro-imaging in particular with data from the Atacama Large Millimeter Array (ALMA). 

O . Methods. We use published radial and proper motion velocities obtained from the shock-excited vibrational emission in the H2 v = 1- 

I S(l) line at 2.122 yt/m obtained with the GriF instrument on the Canada-France-Hawaii Telescope, the Apache Point Observatory, 

O the Anglo- Australian Observatory and the Subaru Telescope. 

^ Results. These data give the 3D velocity of ejecta yielding a 3D reconstruction of the outflows. This allows one to view the material 

C/^ from diff'erent vantage points in space giving considerable insight into the geometry. Our analysis indicates that the ejection occurred 

^ ^^ ^ 720 years ago from a distorted ring-like structure of ~15'' (6000 AU) in diameter centered on the proposed point of close encounter 

of the stars BN, source I and maybe also source n. We propose a simple model involving curvature of shock trajectories in magnetic 

I fields through which the origin of the explosion and the centre defined by extrapolated proper motions of BN, I and n may be brought 

►^ into spatial coincidence. 

^O Key words. Stars: formation - Stars: general - Methods: numerical - ISM: individual objects: OMC-1 

o 

CO 

• 1, Introduction O'Dell 2001; Doi et al. 2002) whereas the second, the so- 

^^ called low- velocity outflow at 1 8+8 km s"\ presents a northeast- 
er The Orion Molecular Cloud 1 (OMC-1) is the nearest and southwest (NE-SW) direction projected onto the plane of the sky 
CN the most studied high mass star-formmg region. The Orion (Genzel & Stutzki 1989; Blake et al. 1996; Stolovy et al. 1998; 
^. Klemmann-Low (Onon-KL) nebula, which harbors an IR lu- Greenhill et al. 1998; Nissen et al. 2007; Plambeck et al. 2009; 
^ minosity of ~10^ Lo (Wynn- Williams et al. 1984), is a region Goddi et al 2009a) 
• '^ of OMC-1 where two spectacular outflows of diflTerent charac- 

/S ter are found emanating from the same confined region of a few The energetic outflows (Allen & Burton 1993; O'Dell 2001) 

^ arcseconds (V = 414+7 AU at the distance of Orion, Menten in OMC-1 have prompted many studies seeking their ages and 

et al. 2007; Kim et al. 2008). One outflow is a wide-angle origins. The age of the fast outflow is reported to be 500 to 

high- velocity outflow, at 30 km s"^ to several hundred km s-\ 600 years (Doi et al. 2002; Bally et al. 2011; Goddi et al. 

oriented in a northwest-southeast (NW-SE) direction (Allen & 2011b). Genzel et al. (1981) reported an age of 3000 years for 

Burton 1993; Chemin & Wright 1996; Schultz et al. 1999; the 18 km s"^ outflow. However, the recent SiO observations 

of Plambeck et al. (2009) carried out with CARMA seem to 

^ „ ~ T ~ u. • J . .u r- J T- TT •• indicate an age «: 3000 years. Nonetheless, our understanding 
Based on observations obtained at the Canada-France-Hawaii ^ , . . ^ , , -^ ^ , ^ , . , - -, , . 
Telescope (CFHT) which is operated by the National Research Council ^^ ^his region has been transformed by high spatial resolution 
of Canada, the Institut National des Sciences de lUnivers of the Centre imaging m molecular hydrogen for which resolution can be bet- 
National de la Recherche Scientifique of France, and the University of ter than 100 mas (Lacombe et al. 2004) m the K-band around 
Hawaii. 2 yum, in the MIR continuum (Shuping et al. 2004) and through 
** Visiting astronomer at the Canada-France-Hawaii Telescope, radio-interferometry of a range of species, most recently methyl 
Mauna Kea, Hawaii. formate, acetone, acetic acid and dimethyl ether (Favre et al. 



H. D. Nissen et al: A 3D view of the outflow in the Orion Molecular Cloud 1 (OMC-1) 



201 la,b; Friedel et al. 2005; Friedel & Snyder 2008; Guelin et al. 
2008), methyl cyanide (Wang et al. 2010) and ammonia (Goddi 
et al. 2011a), including numerous maser observations in SiO, 
OH and H2O (Doeleman et al. 2002; Goddi et al. 2009a,b; Kim 
et al. 2008; Matthews et al. 2010; Wright et al. 1995; Johnston 
et al. 1989; Cohen et al. 2006; Genzel & Downes 1977; Genzel 
et al. 1981; Matveenko et al. 2000; Matveyenko et al. 2007). 
Nevertheless the driving sources of the outflows still remain un- 
certain although it is at least apparent that the high and low ve- 
locity outflows are largely independent of each other and arise 
through difl'erent phenomena. 

While it is evident that convulsive events associated with the 
formation and interaction of embedded high mass stars are the 
origin of the fast outflow in OMC-1 (Nissen et al. 2007), the 
centre of this region presents many candidates (Shuping et al. 
2004; Dougados et al. 1993). A detailed case has recently been 
made that the origin of the high- velocity outflow results from the 
disintegration of a non-hierarchical multiple system (Bally et al. 
2011) where it has been proposed that a close encounter of I, n 
and BN, -500 years in the past, may have triggered the expulsion 
of material (Bally 2008; Zapata et al. 2009; Bally et al. 2011). 
This model is based upon Very Large Array (VLA) proper mo- 
tion measurements of the objects BN, I and n (Gomez et al. 2005, 
2008; Rodriguez et al. 2005). New work of Goddi et al. (201 lb), 
again using the VLA but at the higher frequency of 43 GHz, 
has questioned the involvement of source n while confirming the 
collision at essentially the same position and a similar epoch be- 
tween BN and L These authors propose that source I formed a 
hard binary in the course of the encounter with BN. It has also 
been suggested in a quite diff'erent model that an encounter be- 
tween BN and I, through ejection of BN from the Trapezium 
cluster, may have powered the outflow (Tan 2004). A further 
suggestion is that the nearby radio source SMAl (Beuther & 
Nissen 2008) may be the source of the explosion. In addition, in 
SubmilHmeter Array (SMA) observations of the ^^CO (J = 2-1) 
line Zapata et al. (2011) identified a structure suggested to be 
an expanding bubble with an estimated age of 500-1000 years 
associated with the same central zone. Suffice it to say that the 
nature of the encounter which led to the explosive event 500- 
600 years ago is still a matter of contention but that with a high 
degree of confidence we know that BN, I and perhaps n suff'ered 
a close encounter whose position is known to ~r', irrespective 
of the involvement of source n. 

The low- velocity NE-SW outflow may originate from the 
radio source I (Menten & Reid 1995; Greenhill et al. 2004; 
Beuther et al. 2006). This deeply embedded high-mass young 
stellar object presents a disk (Matthews et al. 2010) perpendic- 
ular to the associated NE-SW ground-state SiO maser emission 
(Goddi et al. 2009a; Plambeck et al. 2009). Alternatively, or in 
addition, source n may play a significant role since it is also 
a massive young stellar object with a possible accretion disk 
(Shuping et al. 2004; Greenhill et al. 2004). More detail on the 
nature and geometry of the low velocity outflow region may be 
found in Nissen et al. (2007), based on K-band H2 observations. 
In Nissen et al. (2007) blue-shifted components are suggested to 
be the NIR counterpart of the radio outflow. Here we concen- 
trate largely upon the fast outflow introducing the low velocity 
outflow separately to form a more complete picture of the dy- 
namics. 

One purpose of the present work is to clarify the geometries 
of the outflows mentioned above through creation of 3D images 
of H2 NIR emission. Such images were first described in Nissen 
(2008). Proper motion measurements (Cunningham 2006; Bally 
et al. 2011) of H2 emission at 2.122 yum with independent data 



for radial velocities (Gustafsson et al. 2003; Nissen et al. 2007) 
are combined here to construct 3D motions as in Nissen (2008). 
The bright emission employed to measure the velocity of the 
disturbance passing through the gas is due very largely to C- 
type magnetic shocks (Colgan et al. 2007; Kristensen et al. 2008, 
2007, 2003; Takami et al. 2002). Shock emission is superposed 
on a nearly ubiquitous background of photon-dominated excita- 
tion, typically at the 10% brightness level. The observations used 
in this work are briefly reviewed in Sect. 2. 

With a model which initially assumes rectilinear ballistic tra- 
jectories, proper motion data (i) enable the determination of the 
approximate location of the expansion centre (already reported 
in Cunningham 2006; Bally et al. 201 1) and (ii) give an estimate 
of the time before present that the explosion took place. The lat- 
ter is described in Sect. 3. In Sect. 4, we construct 3D images 
of the H2 emission in OMC-1. This allows views of the emis- 
sion from any chosen vantage point and also a reconstruction 
of the Orion explosion in 3D. Bally et al. (2011) noted a sys- 
tematic discrepancy between the location of the explosion cen- 
tre and the point (in projection) of close encounter of BN, I and 
n, a discrepancy which in the present work, using a subset of the 
proper motion data and taking into account results in Goddi et al. 
(2011b), appears somewhat heightened. In Sect. 5 we show that 
an inhomogeneous magnetic field could lead to curved trajecto- 
ries of ejecta from a central explosion, removing this important 
deficiency. 

2. Observations: proper motion and radial velocity 
data for molecular hydrogen, H2 

2.1. Proper motion measurements 

Proper motions for the individual features in the OMC-1 out- 
flow were measured using images obtained with the following 
telescopes: i) the Anglo- Australian Telescope (on 13 September 
1992, see Allen & Burton 1993), ii) the Subaru Telescope (on 
the 11 and 13 January 1999, see Kaifu et al. 2000), iii) and 
the Astrophysical Research Consortium 3.5 m telescope at the 
Apache Point Observatory (on 21 November 2004, see details in 
Cunningham 2006; Bally et al. 201 1). In all, 194 proper motions 
were recorded within a 3' x 3' region. Proper motions that are 
used here are in the range of 20 km s"^ up to 120 km s"^ with an 
uncertainty of 10 km s"^. Discussion of errors in proper motion 
data may be found in section 2.3. 

The directions but not the magnitudes of these proper mo- 
tions within the r X r region covered here are shown in Fig. 1. 
On the basis of these proper motion data, Cunningham (2006) 
and Bally et al. (2011) using rectilinear trajectories showed that 
the geometrical centre of expansion of the fast outflow lies at 
0^/2000 = 05^35^^14^.50, ^72000 = -05°22'23r00. In Fig. 1, red ar- 
rows indicate features with position angles which point within 
30° of the centre of expansion. The latter is shown as a deep 
blue dot and is computed from all 194 proper motion data. 



2.2. Radial velocity measurements 

The field of view of radial motions covers the inner 0.1 pc (T 
X r) of OMC-1 and is presented in Fig. 1. Radial motion data, 
ranging up to 50 km s"^ (vlsr), were recorded on December 5* 
2000 at the 3.6 m Canada-France-Hawaii Telescope (CFHT) us- 
ing the Fabry-Perot GriF instrument (Clenet et al. 2002) to per- 
form high resolution IR spectroscopy and obtain radial velocities 
of the H2 emission pixel by pixel from the observed Doppler 
shifts. A detailed description of data and analysis is given in 
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Fig. 1. Field of view covering the inner 0.1 pc (T x 1') of OMC-1 in 
the present work. The grey-scale shows the emission line v=l-0 S(l) 
of H2 at 2.121 //m (observations carried out with the 3.6 m CFHT 
telescope using the GriF instrument, see Gustafsson et al. 2003, 2006; 
Nissen et al. 2007). Green and red arrows, representative of the 71 data 
for which there exist both proper motions and radial velocities, show 
the measured position angles. Those arrows which point with a devi- 
ation of less than 30° from a computed centre of expansion based on 
all 194 proper motion data (Cunningham 2006; Bally et al. 2011) (the 
orange dot - see text) are shown in red and those that show a deviation 
> 30° are in green. Positions of the BN object (0^72000 = 05^35'^ 14U 094, 
Sjiooo = -05°22'22r724), the source n (0^^2000 = 05^35'" 14^.3571, 6j2ooo = 
-05°22'32r719) and the radio source I (0^72000 = 05^^35'" 14^.5 141, 6j2ooo 
= -05°22'3a.'575) (Goddi et al. 2011b) - turquoise dots - are indicated. 
The coordinates of (0,0) are those of the radio source I. Units are arc- 
seconds, where V = 414+7 AU (Menten et al. 2007; Kim et al. 2008). 



Gustafsson et al. (2003, 2006) and Nissen et al. (2007). The spa- 
tial resolution of these data is Qffl^. Radial velocities are shown 
in Fig. 2. 

Seventy-one objects could be unequivocally identified as 
common to the proper motion and radial data sets. Turning to 
the radial data in Fig. 2, the velocities (Isr) cover the range 
-35 km s"^ to -1-50 km s"^ with an average of ~20 km s"^ The 
coordinates of the 88 features observed in the GriF data were 
converted and then compared to the H2 structures in order to 
obtain the same features within the same area. This reduces the 
number of objects to 71. Note that of the 17 features excluded, 
14 were excluded on the basis that the H2 emission was too weak 
to determine an accurate radial velocity and 3 on the basis that 
the objects were not sufficiently resolved that a single radial ve- 
locity could be assigned. Further details of the method used for 
the identification of objects common to both radial and proper 
motion datasets can be found in Nissen (2008). 

2.3. The combined dataset 

It is evident from the position angles shown in Fig. 1 that some 
features cannot reasonably be associated with the outflow. Any 
datum for which the position angle points with a deviation of 
more than 30° from the computed H2 origin, that is, any marked 
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Fig. 2. Radial velocities of 7 1 individual features in the OMC- 1 outflow 
for the same V x V field of view as in Fig. 1. The orientation of the 
image - with the Earth to the top left - illustrates the radial velocities 
associated with each of the objects for which proper motions are shown 
in Fig. 1. Blue and red-shifted emissions are shown as blue and red 
arrows, respectively. The length of the longest arrows corresponds to 
50 km s~^ in the local standard of rest. 



with a green arrow in Fig. 1, was not included in our analysis. 
Thus 14 proper motion measurements were excluded by this cri- 
terion, leaving 57 features in the dataset. We mention in Sect. 
3.2 that this pre-selection has no bearing on our conclusions. 
The choice of 30° arises from a natural dip in the distribution. 
Indeed seven of the discrepant features form a concerted group 
which traverses the field 15'' -20'' northwest of the central re- 
gion. This may indicate that a collimated outflow, unrelated to 
the main OMC-1 flows studied here, crosses the field of view. 

Errors (Icr, see Gustafsson et al. 2003; Cunningham 2006; 
Nissen et al. 2007) in radial velocities and proper motions are 
1 to 2 km s"^ and 10 km s"\ respectively. In this connection 
errors in the proper motions are recorded in Cunningham (2006) 
as 20 km s"^ with a mean value of 15 km s"^ in Bally et al. 
(201 1). In our data for proper motions we use a subset of those 
in Bally et al. (2011) and record an error of 10 km s"^ (Icr). 
This lower value arises because the subset of data which we use 
excludes data (see above) and these exclusions contain a number 
of less accurate values. We also note that sufl&ciently fast shocks 
will dissociate H2. For pre-shock gas of density 10^-10^ cm"^ 
and depending on the local magnetic field strength, this places an 
upper limit of ~70 km s"^ on the maximum velocities of features 
visible in our data (Kristensen et al. 2008). However, some H2 
features with proper motions exceeding this limit are observed 
toward OMC-1. These features are included in the present study. 
It would seem likely that these faster objects impact on material 
which is already in recessional motion, a feature for which there 
is some observational evidence (Nissen et al. 2007). 



3. The OMC-1 outflow : backwards in time 

In this section proper motions of our 57 chosen objects are traced 
backwards in time. Note that here and in all subsequent sections 
up to Sect. 5 rectilinear ballistic motion is assumed. Proper mo- 
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tion measurements record the passage of a disturbance, that is, 
of bullets through the surrounding gas. We consider first if it is 
a good approximation to consider this motion as ballistic. We 
treat the inner bullets shown in Fig. 1 as slower moving versions 
of the fast bullets in the so-called "fingers" (Allen & Burton 
1993) and assume a bullet mass of 5 x lO""^ Mo (Burton 1997; 
Chrysostomou et al. 1997). We use as an example a very well 
characterised shock, 20.7'' west and 6.8'' south of the reference 
star TCC0016, which has been modeled in 3D (Gustafsson et al. 
2010) and imaged at 35 AU resolution (Lacombe et al. 2004). 
The bullet which gives rise to this emission has a velocity of 
50 km s"^ (Gustafsson et al. 2010), typical of bullets in the inner 
region, and thus an energy of -1.25 x 10^^ J. The average emis- 
sion from this region, lozenge shaped and of dimension 0.1" ra- 
dius and 0.5" in length is -10"^ Wm'^sr'^ in the S(l) v=l-0 H2 
rovibrational emission line (see Fig. 11 of Nissen et al. 2007). 
Given these figures, the total energy emitted in the S(l) line is 
~3.3 X 10^^ J if the emission continues for ~700 years as our 
simulations suggest. Emission in the S(l) line makes up -5-10% 
of total H2 emission giving an emitted energy of ~4 x 10^"^ J. 
Allowing for some additional emission lines of other molecules 
this suggests that -4% of the initial energy is lost through radia- 
tion, translating into 2% in the velocity. Thus disturbances in the 
gas involve objects which in general have a considerably greater 
initial kinetic energy than that which is given up in the shocks 
which they generate, justifying the use of a ballistic model as an 
acceptable approximation. 



3.2. The time of the explosion 

Proper motions for the subset of 36 objects in Peaks 1 and 2 were 
run backwards in time, assuming linear ballistic motion. Rather 
than converging on a single spot, these motions converged to a 
ring-like structure shown in Fig. 4 at 720+25 years in the past. 
This time is an upper bound since deceleration has of course 
been ignored in a ballistic model. The computed centre of ex- 
pansion in Bally et al. (201 1) lies at a/2000 = 05^35"" 14^.50, ^/20oo 
= -05°22'23'.'00) and is localized about 7.5" to the north of the 
radio source I (J2000, Plambeck et al. 1995; Goddi et al. 2011b) 
and 6" to the east of the BN object (J2000, Menten & Reid 1995; 
Goddi et al. 2011b). Shown as a deep blue dot in Fig. 4, this is 
found to be located at the centre of the ring illustrated in that 
figure. 

The observed ejecta appear to originate from two caps of 
material, one associated with Peak 1 and the other with Peak 
2, as shown in yellow and green in Fig. 4, respectively, with 
some obvious outliers. There are 5 bullets in the centre, three 
from Peak 1 and two from Peak 2, which as we have seen is 
coincident within measurement errors with the computed centre 
of expansion, based on all 194 proper motion data. If motions of 
the objects are allowed to develop further into the past the bullets 
do not become more confined. Rather, the two caps pass through 
each other and the ring structure is lost within 25 years. These 
characteristics led us to consider the ring structure in Fig. 4 as 
the structure corresponding to the epoch of the explosion. 



3.1. Description of the main OMC-1 outflow 

Constituents of the fast and slow OMC-1 outflows exhibited by 
H2 data within 1' x 1' are shown in Fig. 3. The fast outflow 
is represented by two regions, traditionally called Peaks 1 and 
2 (Beckwith et al. 1978) and labelled as such in Fig. 3. The 
slow outflow, first detected in H2 in Gustafsson et al. (2003), 
and which we call Region B (see Fig. 3), is composed entirely 
of blue- shifted shock zones (Nissen et al. 2007). 
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Fig.3. The naming conventions for Peak 1, Peak 2 (Beckwith et al. 
1978) and Region B towards OMC-1 superposed on the image of the H2 
v=l-0 S(l) line at 2.121 yt/m. The coordinates of (0,0) are those of the 
star TCC0016 (aj2ooo = 05^35'" 14!909, 6j2ooo = -05°22'39'.'31), Vannier 
et al. (2001). The red dot indicates the position of source L Units are arc- 
seconds, where V = 414+7 AU (Menten et al. 2007; Kim et al. 2008). 
Figure taken from Nissen et al. (2007). 




Fig. 4. Disposition of the Peak 1 and Peak 2 objects in yellow and green, 
respectively at the epoch time zero of the explosion seen in the line-of- 
sight (i.e. perpendicular to the plane of the sky xy). Thirty-six features 
are displayed in this figure. The dark dot is the computed expansion cen- 
tre (Nissen 2008; Bally et al. 201 1). The ring-like structure which forms 
when proper motions are run backwards in time is highlighted in violet. 
This structure is the starting point for the 3D-simulations illustrated in 
subsequent figures. Units are in arcseconds. 



The observed ring- structure has a diameter of about 6000 AU 
(~ 15'0- The empty lane between the centre and rim seen in Fig. 4 
shows that the structure is not a sphere given that we are observ- 
ing clumps of dense shocked gas. The structure is robust to in- 
clusion of errors and may be clearly discerned with use of all 194 
proper motion data without introduction of any of the selection 
criteria adopted here. Results also show that the ring structure is 
largely populated from objects within ~0.1 pc of the expansion 
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centre. Thus we do not propose this structure as the precursor of 
the very fast moving bullets which make up the outer fingers of 
emission (Allen & Burton 1993). Figure 2 of Bally et al. (2011) 
illustrates these points. 

Further in connection with the ring, we note that the ring 
is well described as circular and this suggests it is not tilted in 
the plane of the sky. Moreover it forms by a remarkable coin- 
cidence of individual and numerically very disparate values of 
proper and radial motions. An apparent lack of structure up to 
-650 years in the past leads to a rather abrupt appearance of 
the ring structure (± <25-50 years) running further backwards 
in time. Beyond 720 years in the past the ring structure rapidly 
disappears. 

The apparent origin of so many bullets in a well defined 
structure near the computed centre at one time lends consider- 
able weight to a single explosion hypothesis to describe the in- 
ner part of the fast outflow lobes. However there are aspects of 
the data that highlight the limitation of our present simple inter- 
pretation. For example, three Peak 2 bullets in Fig. 4 appear to 
originate from Peak 1 . In addition, as noted above, the explosion 
centre in Figs. 1 and 4 diflfers from the position of closest en- 
counter of I, BN (and n) (Bally et al. 2011; Goddi et al. 2011b; 
Zapata et al. 2011), a point which we discuss in more detail in 
Sect. 5. 

At all events our analysis indicates that the ejection occurred 
< 720 years ago, which is in reasonable agreement with the 
range of dates of 500-560 years deduced in previous studies 
(Bally et al. 2011; Goddi et al. 2011b; Gomez et al. 2008, 2005; 
Rodriguez et al. 2005; Zapata et al. 2009) and data show, as in 
Bally et al. (2011), that the close encounter of BN, I and n is a 
feasible explanation for the Orion explosion. 

4. A 3D view of the OMC-1 outflow 

Our purpose here is to combine proper motion and radial data 
to create a 3D model of the Orion explosion as observed in 
H2 emission. We remark first that proper motion measurements 
record the passage of a shock through the surrounding gas, 
where this disturbance, marked by H2 excitation, is generated 
through material - bullets - flung out at high velocity due to an 
explosive event. Radial velocities obtained from GriF data how- 
ever record the actual motion of the H2 itself and not directly of 
the bullets. We base all 3D constructions that follow on the as- 
sumption that the velocity in the radial direction can be treated as 
a measure of the radial velocity of the disturbance and thus that 
the gas that is emitting has been accelerated to the full shock ve- 
locity. In fact C-type shock models show that the peak emission 
in the gas is reached when the gas is moving typically at a ve- 
locity which is 5 km s"^ or more below the velocity of the shock 
and thus the radial motion of the hot H2 is not a precise measure 
of the radial component of the motion of a bullet. This velocity 
diflTerence arises because it is the conversion of kinetic energy 
into internal energy which generates the hot H2. However since 
motion is concentrated in the plane of the sky the associated sys- 
tematic error in the estimate of 3D velocities is <2-3 km s"^ in 
values which may be typically 50-60 km s"^ (Kristensen et al. 
2008). This small systematic eflTect has been ignored here. 

4.1. The fast outflow: Peaks 1 and 2 

Based on the above conclusion that the ejecta originate from a 
ring-like structure, and not some more spherical distribution, the 
structure shown in Fig. 4 is taken to mark the epoch of the ex- 
plosion and all bullets that form Peaks 1 and 2 are assumed to lie 



in the plane of the sky at this epoch. Motion now runs forward 
in time, including both radial and proper velocities, evolving to- 
wards the present and tracing out the formation of the lobes of 
the explosion. The 3D structure of the H2 bullets 720 years later 
is shown in Fig. 5. Snapshots at 40 year intervals, including slow 
outflow objects, are shown in Appendix A (Fig. A.l)^ 




Fig. 5. Disposition of Peak 1 and Peak 2 objects in yellow and green 
respectively at the present epoch seen in the line-of- sight. Thirty-six 
features are displayed in this figure. This image is a reconstruction from 
Fig. 4 after 720 years. The dark blue dot is the computed expansion 
centre, as earlier. Larger and brighter objects are closer. Units are in 
arcseconds. 



The appearance of Peak 1 and 2 bullets at a viewing angle of 
45° from below the plane of the Galaxy is illustrated in Fig. 6. 
This figure shows that Peaks 1 and 2 regions are each bent away 
from the Earth at an angle of ~ 15° out of the plane of the sky. The 

two lobes can be seen to lie at an estimated value of ~150°^^^o 

-20 

rather than opposed at 180°. The true angle may be smaller than 
150° because of the observational selection of proper motions 
and radial velocities. Since both flows are angled away from the 
Earth, the efl'ect does not arise through dust obscuration. Note 
that the eff'ect is already apparent in the radial data shown in 
Fig. 2. The presence of an inhomogeneous magnetic field could 
be responsible for this phenomenon, as discussed in Sect. 5. 

4.2. The slow outflow: Region B 

Region B is a continuous or punctuated outflow running from 
the northeast (NE) to southwest (SW), located southwest of ra- 
dio sources I and n (see section 3.1 and Fig. 3). The reason for 
the choice of a continuous outflow for Region B rather than an 
origin in a single explosive event is given below: at all events the 
nature of the Region B outflow is quite distinct from the explo- 
sive fast outflow to which it lies roughly at right-angles in the 
plane of the sky. There is clear observational evidence, set out in 
Nissen et al. (2007), that the Region B outflow originates from 
a central zone close to or coincident with source I. However it 
remains uncertain whether the origin of the outflow is source 
I, source n or some other unknown object. Thus, Nissen et al. 
(2007) have suggested that Region B is the IR counterpart of 
a continuous outflow detected in the radio from either source I 



^ Animations allowing a visualization from this and other vantage 
points can be obtained upon request to the authors. 
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Fig. 6. Disposition of Peak 1 and Peak 2 objects in yellow and green re- 
spectively at the present epoch viewed from an angle of 45° beneath the 
Galactic plane, that is, the xz plane. Earth lies in a direction vertically 
upwards. The plane of the sky is shown corresponding to z = 0, passing 
through the computed expansion centre. An angle of ~ 150° between the 
axes of the Peak 1 and 2 flows may be seen. Units are in arcseconds. 



or n (Greenhill et al. 1998, 2004; Matthews et al. 2010; Goddi 
et al. 2009a; Plambeck et al. 2009; Wright et al. 1996; Genzel 
et al. 1981). Recent v=0 SiO maser data (Goddi et al. 2009a; 
Plambeck et al. 2009) show very clearly the association of this 
extended SiO emission with source I. The northeast-southwest 
orientation of the SiO emission provides further evidence that 
the Region B outflow is associated with source I on the geomet- 
rical grounds of the flow orientation. In conclusion, since there 
is good evidence in Nissen et al. (2007) and from more recent 
SiO masers observations, we feel justified in treating Region B 
as a counterpart of a continuous (or punctuated) outflow arising 
from the central zone close to the origin or at the origin of the 
explosion which gave rise to Peaks 1 and 2. 

The difl&culty accompanying this characterization of Region 
B is that the BN, I (n) encounter took place -700 years ago, us- 
ing the figure derived from our data. The accompanying violence 
which gave rise to Peaks 1 and 2 at that time would disrupt any 
outflow from I itself or from the near vicinity of I, if only since 
outflow B is of much lower total energy content. Thus we require 
that the outflow B is no older than -700 years and it would seem 
necessary to put on one side the age estimate of Genzel et al. 
(1981) of 3000 years and rather adopt the value of a few hun- 
dred years consistent for example with data in Plambeck et al. 
(2009). This is itself consistent with our finding that extrapola- 
tion of ballistic motion of objects in Region B backwards in time 
for the full 720 years (see section 4.1) results in some objects in 
the Region B outflow passing through the central zone and re- 
treating out the other side. Our inference is that in these cases 
the H2 emission traces more recently emitted objects in a con- 
tinuous outflow from the central zone. Using this model, data 
for Region B are now combined with those for the fast outflows 
already illustrated. 

4.3. Peaks 1, 2 and Region B objects 

We now seek to combine all objects in our sample. We first ex- 
plain how Region B objects can be incorporated, given the con- 
siderations in the preceding section. As noted, the Region B out- 
flow is of quite diflferent character to that of the Peaks 1 and 



2. The latter originate from a single explosive event whereas 
Region B objects constitute members of a much more familiar 
form of continuous or punctuated outflow associated with star 
formation. The Region B outflow is therefore treated as emanat- 
ing from a single zone with continuous replenishment as time 
passes. If we take the present positions of Region B objects and 
reverse them in time, such a reversal would therefore endure for 
a length of time which is diflTerent for every object. This set of 
times is given by the times required for objects to find their way 
back to the zone of origin of the outflow assuming rectilinear, 
ballistic motion. In our simulation the zone of origin is identified 
as any point on a line bisecting Peaks 1 and 2. Thus if we now 
move forwards in time, tracing the spatiotemporal evolution of 
the Region B outflow, new objects are created at this set of times, 
lying between 720 years in the past and up to close to the present. 
Region B objects as a whole involve objects with radial veloci- 
ties in the range -24 km s "^ to -8 km s "^ and proper motions 
from 28 km s "^ to 102 km s "^ There is a notable diflTerence be- 
tween these velocities and values for 22 GHz H2O masers which 
are clustered around 18 km s"^ (Genzel et al. 1981). H2O masers 
form where the temperature is very much lower and the density 
very much greater than in the region where H2 emission is ob- 
served (Yates et al. 1997). Hence, we expect H2O masers and H2 
features to trace diflTerent portions of the flow. 

The present day appearance of the objects associated with 
Peak 1, Peak 2 and Region B are shown in Fig. 7 from diflTer- 
ent vantage points: a) as seen from Earth, b) as seen from the 
direction of the Sagittarius arm of the Galaxy^ and c) as seen 
from beneath the Galactic plane. More specifically. Fig. 7c and 
to some extent 7a, show how Region B protrudes towards the 
Earth (along the negative z axis) and points towards the west. 
Fig. 7b also shows that Region B is angled strongly south, not- 
ing that this eflTect may be exaggerated through the omission of 
blue-shifted objects with proper motions less than 10 km s"^ 

Thus Region B represents a flow in a south-westerly direc- 
tion which is angled to some small extent towards the Earth. 
This geometry can be rationalized if source I is the source of the 
Region B outflow, as posited in Nissen et al. (2007). Consider 
the plane defined by the direction of proper motion of source I 
and the normal to the sky. According to results in Goddi et al. 
(201 lb) this is a plane tilted 20° south to the horizontal with the 
motion of I pointing south-east. Source I is recorded in Goddi 
et al. (2011b) as moving at 11.5 km s"^ which is reported to be 
the full velocity of source I with respect to the nebula. Source 
I is thus moving in the plane of the sky. An outflow from I 
will therefore appear to be projected in a manner which com- 
bines the motion of I and the radial outflow velocity of Region 
B, whose average is 18 km s"\ a figure about which values are 
quite tightly clustered (see Fig. 5 of Nissen et al. 2007). This will 
bend the flow to the south and make it a little less westerly com- 
pared with the orientation of the outflow in the moving frame of 
source I. The same consideration would apply to outflows asso- 
ciated with extended SiO maser emission in v=0 (Goddi et al. 
2009a; Plambeck et al. 2009; Goddi et al. 2009b; Cho et al. 
2005) as alluded to in Plambeck et al. (2009). In this connection 
the orientation of the Region B outflow, as seen most clearly in 
Fig. 7c, corresponds well with the axis of the SiO v=0 outflow 
recorded in Goddi et al. (2009a) as mentioned in section 4.2. In 
both cases these show a position angle of -230° measured north 
to east from the axis defined by the fast outflow. The direction 
of the Region B outflow would also appear to be at right-angles 



^ Orion lies in the Galactic plane and the Sagittarius arm is at right 
angles to the Earth - OMC-1 direction. 
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Fig. 7. Peak 1, Peak 2 and Region B objects at the present day in yellow, green and turquoise, respectively. The deep blue dot is the computed 
expansion centre. Units are in arcseconds. a) The view as seen in the line-of-sight, b) the view as seen from the Sagittarius arm of the Galaxy 
(normal to the yz-plane), c) the view as seen from beneath the Galactic plane (xz). 



to the orientation of an edge-on disk around I, whose presence is 
inferred from SiO v=l,2 maser observations in Matthews et al. 
(2010) as mentioned in the introduction. The Region B outflow 
will be sprayed out to some extent due to dispersion about values 
of the radial velocity and proper motion velocities. Objects will 
therefore be disposed relative to the source of the outflow de- 
pending both on the time of ejection of the bullets and their indi- 
vidual proper and radial motions. The protrusion of any Region 
B bullet out of the plane of the sky and towards the Earth will be 
tan" ^(radial velocity/proper motion velocity). This is -20° using 
a typical value of radial velocity of 18 km s"^ and proper motion 
velocity of 50 km s"^ 



Above for simplicity we have arbitrarily used the expansion 
centre for the origin of the Region B outflow. If however we were 
to use source I, say, then features observed in Region B would 
have to be introduced slightly later into the model. This is due 
to the geometrical eflTect of the shift of the SE-NW line bisecting 
Peaks 1 and Peaks 2 that now runs through source I instead of 
the centre of expansion. Accompanying changes in the positions 
of Region B objects in the 3D images do not exceed 3'\ 



5. Magnetic deflection of the ejecta in the OMC-1 
explosion 

A detailed and convincing model has been developed (Bally 
et al. 2011; Goddi et al. 2011b; Gomez et al. 2008, 2005; 
Rodriguez et al. 2005) in which the fast outflow in Orion was 
caused by the close encounter of sources BN, I and n -500 years 
in the past, notwithstanding the doubts recently raised about the 
involvement of source n. A potential problem with the hypoth- 
esis, noted in Bally et al. (2011) and to which we have already 
referred, is that the proper motions of the 194 ejecta recorded 
in that paper, produced as rectilinear, intersect at a centre which 
is displaced from the known position of the explosion by ~5'' 
to the northeast. Since the point of intersection has itself an un- 
certainty of 5'\ Bally et al. (2011) concluded that the position 
of the explosion and of the intersection were not difl'erent at a 
level to cause too much concern. However this remains a nag- 
ging problem and it is heightened in the present work where the 
centre of the ring in Fig. 4 is 5.5'' south and 2.5'' west, that is, 
6" displaced from the BN, I, n centre of encounter or 7" dis- 
placed relative to the centre presented in Goddi et al. (2011b). 
In addition, the present apparent ejecta centre in Fig. 4 is rather 
closely spatially confined by 5 bullets occupying a region of di- 
ameter <3" with the BN, I, (n) centre of encounter lying on or 
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very near the rim of the circular structure in Fig. 4. Errors in the 
absolute position of the BN, I, (n) encounter are small by com- 
parison (~r', Goddi et al. 2011b) and do little to alleviate this 
problem. 

A further issue is that the outflows are unexpectedly bent 
away from us in the plane of the sky in both Peaks 1 and 2, as 
recorded in section 4.1 and figure 6. Such a geometry does not 
arise naturally from the extensive models in Bally et al. (2011) 
or Goddi et al. (2011b). Both these difficulties may be simulta- 
neously resolved if the ejecta have been moving since the time 
of the explosion in an inhomogeneous magnetic field showing an 
organised structure on the scale of the inner 0.1 pc of the outflow. 
As we have noted, shocks and the accompanying vibrationally 
excited H2 are generated through the rapid passage through the 
medium of masses of material - the bullets. The physical nature 
of these bullets is unknown but it is reasonable to assume that 
by whatever of the three explosion mechanisms proposed for ex- 
ample in Bally et al. (201 1), the material of the bullets may have 
achieved substantial ionization. The trajectories of these ejecta 
will then become curved in the presence of a magnetic field gra- 
dient. At this juncture we should add that it is possible that the 
deflection is purely hydrodynamic. Specifically, the dense OMC- 
1 dust ridge runs mainly north- south. Thus outflow components 
moving north-west may run into a density gradient which natu- 
rally will tend to deflect the ejecta towards the west. We do not 
pursue this line of argument further, if only through the lack of 
observational constraints. 

To return to the magnetic model, deflection of a plasma jet in 
a magnetic field gradient is a well-known phenomenon in plasma 
physics and engineering (e.g. Demichev et al. 1965; Timoshenko 
2008). The phenomenon does not appear to have been much con- 
sidered in the astrophysical context in contrast to magnetic col- 
limation or magnetic centrifugal acceleration (e.g. Bogovalov & 
Tsinganos 1999; Banerjee & Pudritz 2006). The most simple 
manner to appreciate the phenomenon is through the magnetic 
pressure, BV2yUo (Demichev et al. 1965). If this increases radi- 
ally, in the direction of the radius vector r, then the excess pres- 
sure on the opposite or outer side of any one of the ejecta will 
be greater than on the inner side, giving rise to motion in a cir- 
cular path. This approach glosses over the physics involving the 
detailed motions of the electrons and ions but allows an estimate 
to be made of the radius of curvature, r, of the path of plasma 
ejecta. The relatively small eff'ect of deflection in the direction 
of r X B, that is, perpendicular to the radius vector, is ignored 
here (Timoshenko 2008). 

For simplicity, imagine a magnetic field whose structure is 
a continuous small segment of the surface of a sphere. The true 
point of ejection, which we take to be the BN, I (n) centre of 
encounter, lies at some position within the sphere defining the 
segment. The segment is taken to be convex in our direction with 
the field increasing with radius. As a mass of ejected material en- 
counters this field the diff'erential magnetic pressure across any 
one mass of material will cause its path to follow a continuous 
geodesic, moving either northwest (Peak 1) or southeast (Peak 
2) showing curvature both in the plane-of-the-sky and in a di- 
rection away from the observer, that is, the line-of-sight. Both 
motions will have the same radius of curvature. If we take tan- 
gents to these paths, equivalent to taking rectilinear trajectories, 
these tangents will meet at an apparent origin north of the true 
origin of the explosion projected onto the plane-of-the-sky. The 
present positions will also be bent away from the observer. 

We now seek to show that the displacement of the apparent 
origin from the BN, I (n) centre of encounter can be 6'' - 7'' given 
the known or likely properties of the system, and further that the 



angle at which the outflows are bent away from the observer is 
-15° under the same set of physical conditions. We derive first 
some simple equations which govern the motion of the ejecta. 
If the diff'erence between the magnetic field on the outer side of 
one of the ejecta and on the inner is AB, with B outer > B inner and 
the material eff'ectively presents a surface of radius rb, then the 
net force sustaining the circular geodesic path of radius r will be 
7rr^AB^/2yUo where 7rr^AB^/2yUo = mbV^/r where mb is the mass 
of any bullet. Thus, 



2yUombV^ 

AB^Trr^ ' 



(1) 



The form of Eq. 1 is essentially the same as that derived in 
Demichev et al. (1965) given that the gradient of the field is lin- 
ear within the mass of material. 

If the path sweeps out an arc covering an angle 6 then the 
mass will deviate from an extrapolated position based on a 
straight line trajectory to its true position by a distance, s, re- 
lated to 6 through cos(6) = (1-s/r). Thus for small s/r we find 
6 ~ (2s/r)^/^. The distance travelled is Or and therefore the time 
travelled is 6r/\ given ballistic motion. Hence the time of travel, 
t, is given by 



t = 



2(yUombS/7r) 
ABrb 



1/2 



(2) 



A cartoon illustrating this model is shown in Fig. 8. 

If we make the simplifying assumption that projectiles tend 
to be clustered around values of similar dimension and density it 
is evident that the displacement, s, will be similar for all ejecta. 
It is clear from Fig. 4 that there are some outliers and this can 
be attributed to a breakdown of this assumption, if not to other 
over-simplified aspects of the model such as a single value of net 
magnetic field. 

The time t is taken to be 720 years from the present work. 
This may be compared with the figures of 500 and 560 years 
reported in Bally et al. (2011) and Goddi et al. (2011b) respec- 
tively. We take the typical radius of a bullet to be 1.0'', esti- 
mated from the size of the perturbation created through the pas- 
sage of the shock and the mass mb to be a 5x10""^ Mo (Burton 
1997; Chrysostomou et al. 1997). Note that values of masses and 
velocities are consistent with the total energy of the outflow of 
-10^0 J (e.g Bally et al. 2011; Chrysostomou et al. 1997). These 
parameters imply a density of nn + 2nH2 ~ 6x10^ cm"^, taking 
account of He, comparable with estimates that may be found for 
example in Tang et al. (2010) and consistent according to those 
authors with a magnetic field of >9 mG. The observed displace- 
ment between the apparent centre and the actual centre of explo- 
sion is -4x10^"^ m (6.5''). From Eq. 2, it follows that AB would 
have to assume the value of ~6 mG to achieve this displacement. 

As mentioned above ejecta will also be bent away from the 
observer by the magnetic pressure engendered by AB. Simple 
geometrical considerations, given motion on a sphere, show that 
the angle at which ejecta are deflected is given by tan"^(s/d) 
where d is the distance that ejecta would have traveled from the 
apparent explosion centre in the absence of any magnetic deflec- 
tion. Using values from above, it follows that this angle equals 
14°. This implies an angle of 152° between the flows, agreeing 
very well with the albeit approximate value of 150° shown in 
Fig. 6. Thus the magnetic deflection model has the advantage 
of ofl'ering an explanation for both the problem of the apparent 
lack of coincidence of the origin of the outflowing bullets with 
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Fig. 8. Cartoon illustrating the model adopted for the deflection in a magnetic field for an individual bullet of mass mb and radius rb. The bullets 
which create the shocks are assumed to experience an excess magnetic pressure such that ^outer > dinner (top left). The bullets move with a velocity 
Vb along a curved path due to the excess magnetic pressure which they thereby feel (bottom left). The figure bottom right shows the displacement 
in the plane of the sky, s, of a bullet between that arising from a rectilinear motion in the absence of a net magnetic pressure and the curved path 
in the presence of a magnetic pressure. The angle Q, see the derivation of Eq. 2, is also illustrated. 



the BN, I (n) encounter centre and also the deflection backwards 
out of the plane of the sky of Peaks 1 and 2. 

Data in Houde et al. (2004); Vaillancourt et al. (2008); Tang 
et al. (2010) demonstrate the presence of a highly ordered mag- 
netic field covering the region of interest in Fig. 1 at a posi- 
tion angle -130° east of north. With regard to absolute values 
of the field, Tang et al. (2010) argue that a (tangential) mag- 
netic field of anything between 3 and 30 mG may be required 
to confine ejecta. These considerations suggest that an ordered 
field of several mG is appropriate for the gas through which 
the ejecta travel. However within the cavity created by expul- 
sion of material due to the explosion, the magnetic field would 
be likely to be very low since the field would be expelled with 
the ionized material. Note that a marked absence of material is 
observed in the region around the ejection centre. This is illus- 
trated for example by the absence of dust emission around the 
explosion centre in continuum maps in Tang et al. (2010). It is 
also illustrated by the absence of H2 emission and warm dust 
emission in images in Lacombe et al. (2004) where H2 emis- 
sion should otherwise be generated through the high UV flux 
from the Trapezium in this region were there material present at 
densities of a few times 10^ cm"^ or higher. Thus we propose 
a fragment of a spherical cavity with material outside the cav- 
ity maintaining a field of several mG and a small field within 
the cavity. This would give rise to the magnetic gradient which 
is a necessary ingredient of the above model for deflection of 
the trajectories of ejecta. In this connection magnetic fields de- 
rived from OH maser Zeeman splittings vary very considerably 
- between -3.5 mG and -1-16.3 mG - within the central arcsecond 



around the position of radio source I and do not show a clear pat- 
tern (Cohen et al. 2006). These fields are local to each maser spot 
and may be wound up by local turbulent eddies (Field 1982). 

Notwithstanding the broad success of this magnetic deflec- 
tion model we have omitted to include magnetic deflection in 
tracing out the motions of the bullets in any of the relevant fig- 
ures in this paper. This choice has been made in order to maintain 
the results as strongly observational based as feasible. The intro- 
duction of magnetic deflection would necessitate assumptions 
about the mass, magnetic field and bullet density distributions. 

6. Conclusions 

Combining proper motion and radial velocity data it has proved 
possible to create 3D images of the structure of OMC-1 in vi- 
brationally excited H2 in the inner 0.1 pc. These graphic images 
are based directly on observational data and require few assump- 
tions for their creation. The most significant assumptions are that 
the H2 bullets lie in the plane of the sky at the epoch of the ex- 
plosion and that the motion is ballistic. This latter assumption 
may readily be relaxed without any qualitative implications for 
our understanding. 

The magnetic deflection model introduced here, simple as it 
is, removes any lingering doubt that the outflow was in some way 
triggered by the conjunction of BN, I (n) at an epoch between 
500 and 700 years in the past. This model is able to reconcile 
the systematic departure of the apparent explosion centre, based 
on linear extrapolation of proper motions, with the position of 
the BN, I (n) encounter, using parameters which appear quite 
appropriate to the region of interest. 
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Results presented here suggest that the explosion impacted 
on a pair of distorted arcs of molecular material forming a ring 
~15'' (6000 AU) in diameter. We also see how Peaks 1 and 2 
are swept away from the observer into the plane of the sky by 
-15°. Images further reveal the disposition of the slow blue- 
shifted outflow, Region B, which is found to be angled to the 
south and west in a manner which does not appear to relate to 
the fast outflow but may be connected with the motion of I if it 
is the source of the slow outflow. 

In conclusion, while discrepancies remain with regard to the 
details of the close encounter of massive bodies within the core 
of OMC-1 and the overriding contributions to the mechanism 
which led to the explosion, a remarkably consistent phenomeno- 
logical picture now emerges at least with regard to the fast out- 
flow forming Peaks 1 and 2. 
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Appendix A: 3D view of the OMC-1 outflow 

The following figures allow the 3D reconstruction of the OMC-1 
outflow observed in the H2S(l)v=l-0 line from the expansion 
centre < 720 years ago to the present time. The material may be 
viewed from diflTerent vantage points in space. Complete anima- 
tions may been obtained upon request to the authors (D. R, H. N. 
or C. E). 
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Fig. A.l. 3D view of the OMC-1 outflow from the expansion centre < 720 years ago to the present time. The blue dot is the computed expansion 
centre (CK72000 = 05^^35'" 14^.50, (^72000 = -05°22'23'/00). Yellow, green and turquoise points show the disposition of the Peak 1, Peak 2 and Region 
B objects through time. The snapshots are at intervals of 40 years, xy is the plane of the sky. Notwithstanding the magnetic deflection model 
introduced in section 5, we have for simplicity and to avoid model dependence used rectilinear trajectories in this and all other relevant figures. 



